Fitting CoGeNT Modulation with an Inelastic, Isospin- Violating Z' Model 
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We reanalyze the annual modulation data observed by CoGeNT experiment and show that if 
the annually modulated anomaly detected by CoGeNT is induced by collision between dark matter 
particle and nucleus, it can be fitted by a Z' model with inelastic dark matter and isospin-violating 
interaction, and the constraint from XENONIOO can be avoided. This Z' model is strongly con- 
strained by collider physics that the upper bound of the mass of Z' is around twice of the mass of 
dark matter. 

PACS numbers: 



I. INTRODUCTION 

A lot of astronomical evidence have already shown that 
more than one fifth of the universe is composed by mat- 
ter which is not baryonic and does not significantly in- 
teract with light, the nature of dark matter (DM), how- 
ever, remains mysterious. The most compelling version of 
DM is that it is composed by weakly interacting particles 
(WIMPs) , which offers to understand the relic abundance 
of DM as natural consequences of the thermal evolution 
of the universe. 

A number of experiments are searching for a signal 
from WIMPs scattering on nuclei with underground de- 
tectors. Among them, CoGeNT collaboration, employing 
P-type point contact germanium, is able to detect signals 
with sub-keV nuclear recoils due to their low electronic 
noise. In early 2010, around 100 events above the ex- 
pected background were reported by CoGeNT with ion- 
ization energy in the range of 0.4 to 1.0 keVee (electron 
recoil equivalent energy) [l| with a 56-day exposure in 
330 grams germanium inside the fiducial volume. These 
events can be interpreted as elastic scattering of a 5 — 10 
GeV WIMP on nuclei with ~ lO"''" cm^ cross section. 
Most recently, the CoGeNT collaboration published their 
new results of 442 live days data ^ , which confirmed the 
event excess observed in their previous report and showed 
a 2.8(7 indication of an annual modulation of the event 
rate. However, as noticed in the literature [sl, the best fit 
point of the modulation signal is excluded by the unmod- 
ulated rate. On the other hand, the XENONIOO collabo- 
ration alsopublished their new results from a 4843 kg day 
exposure [j], obtaining three events passing all the cuts 
with an expected background of 1.8±0.6 events. The 
large exposure and fiducial volume render XENONIOO 
the most stringent constraint for heavy WIMPs ('-^lOO 
GeV). However, the specific energy resolution at low en- 
ergy region gives XENONIOO a strong rejecting power for 
light WIMPs (~10 GeV) as well [ll[, which appears to 
rule out the light WIMP interpretation of CoGeNT data. 
However, there are several ways that this conflict can be 
alleviated. From experimental side, there might be issues 
with scintillation efficiency of liquid xenon or the quench- 



ing factor of the germanium in the region of interest . 
From the theoretical point of view, isospin violatin g iii - 
teraction between WIMP and nucleons is proposed |15| . 
where, in particular, if An/Ap « —0.7, the constraint 
from xenon-based experiment on germanium based ones 
would be strongly weakened by a factor of around 20, 
where Ap and An are the scattering amplitude of WIMP 
on proton and neutron, respectively. 

In summary, it faces two major challenges to inter- 
pret the result of CoGeNT as WIMP signals. One is 
the tension between the unmodulated rate and the mod- 
ulation from itself, and the other is the conflict from 
XENONIOO. The bound on WIMP-nucleon cross section 
of XENONIOO reduces sharply with the mass of WIMP 
due to that the nuclear recoil energy is proportional M|, 
in the case of low mass dark matter, where Mr, is the 
mass of DM, and XENONIOO is not very sensitive in low 
energy region. For example, the upper bound on scat- 
tering cross section between WIMP and nucleon weakens 
by a factor of 20 from Md = 7 GeV to Md = 6 GeV. 

In most prevailing models, the DM candidate has only 
elastic scattering with nuclei, such as the LSP DM in su- 
persymmetry models [5i] or the KK DM in models with 
extra dimensions 0. However, motivated by the annual 
modulation signal from DAM A/LIBRA experiment f?], 
more attentions have been drawn to the idea of inelas- 
tic DM models [8|. After the modulation signal of Co- 
GeNT was announced, people started to think about in- 
terpreting the CoGeNT result as signals induced by in- 
elastic WIMP 11, 111 . In this work, we try to fit the Co- 
GeNT modulation data using a leptophobic Z' portal 
with isospin-violating couplings to quarks in the Stan- 
dard Model (SM). Furthermore, the DM in this model 
is a real scatter with a keV energy gap from its excited 
state, which is generated by a soft breaking of a global 
U(l) symmetry. Apart from DM direct detection exper- 
iments, this model is also strongly constrained from the 
process of monojet plus missing transverse energy (MET) 
in Tevatron and LHC. 

The rest of this work is organized as the following. In 
Sec. II, we introduce the basic formulas of DM direct 
detection for future use. In Sec. HI, we derive a model 
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independent formula to translate the differential rate of 
WIMP scattering from one experiment to another in the 
region of low nuclear recoil energy and use it to motivate 
the idea of inelastic WIMP. In Sec. IV, we present ba- 
sic properties of the model and use the model to fit the 
modulation detected by CoGeNT. In Sec. V, we discuss 
the collider constraint on this model. We summarize in 
Sec. VI. 



II. ELASTIC AND INELASTIC WIMPS 

The recoil energy Er of the nucleus hit by WIMP with 
velocity v in lab frame can be written as 
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where 9 is the scattering angle in the center-of-mass 
frame of the WIMP and nucleus and S is the energy 
gap between WIMP and its excited state, ^da is the 
reduced mass of WIMP and nucleus. Since the mass of 
WIMP we are interested in in this work is around 10 
GeV, the values of (J-da in Germanium and Xenon are 
similar to each other. Therefore, from Eq. ([T|), one can 
see that the dependence on target nucleus relies on niA, 
and the dependence is quite simple so that one can trans- 
late the differential energy spectrum from one detector to 
the other. 

The differential scattering rate of WIMP with mass 
Mo and velocity distribution /(v) can be written as 
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where Ma and Ej. are the mass and the recoil energy of 
the target nucleus, and Wmin is the minimal velocity for 
the WIMP to generate certain nuclear recoil energy Er, 
which can be written as 
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where 5 is the energy gap between DM and its lowest 
excited state, (5 = in the case of elastic scattering. 

The differential cross section of WIMP and nucleus in 
Eq. ^ can be written as 
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where hdn '^s the reduced mass of WIMP and nucleon. 
On = {<Jn + crp)/2 is the average scattering cross section 
between WIMP and nucleon and 
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2 G isotopes 



2r, [Z cos 6'jv + (A, - Z) sin 6ljv]^ , (5) 



where are relative abundances of isotopes, tan^jv = 
An/Ap. In Eq. Q, F{c^) is the form factor of nucleus. 



and in this work, we assume Fn = Fp = F and use form 
factor proposed by Helm [l^l that 



(6) 



where ji is the first order spherical Bessel function and 
r„ is an effective nuclear radius can be written as 
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where c ~ {1.2iA^/^ - 0.60) fm and s = 0.9 fm. 

The differential scattering rate in Eq. ([2]) depends 
strongly on the velocity distribution of dark matter halo. 
In this work we assume a standard halo model with the 
isotropic Maxwellian velocity distribution: 



cxp[ 
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where ve is the velocity of the earth, Vcsc is the local 
Galactic escape velocity and k is the normalization factor 
which is 



2\3/2 



erf 



1-1/2 



Vo 



(9) 



where erf is the error function and Ucsc is the escape ve- 
locity of DM in local galaxy which is about 544 km/s 

The velocity distribution of DM is formed due to grav- 
itational effect. If the collision between DM and other 
regular matter in the universe is dominated by inelastic 
scattering, the kinetic energy of the system gets smaller 
during the collision. Therefore, if there are two compo- 
nents of DM in the same halo, one is elastic and the other 
is inelastic, one may expect that the average velocity of 
the inelastic DM is smaller than the elastic one. There- 
fore, in this work, we choose a relatively small vo which 
is 200 km/s. 

/(ve) is modulated annually due to the rotation of 
the earth around the sun, which then induces the annual 
modulation of the scattering rate between WIMP and 
nucleus. The velocity of the earth can be parameterized 
as the following flill: 
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(10) 



where v^^^ — 232 km/s and Wg"* = 15 km/s. 

Due to the sensitivity limit of the detector and the non- 
linear velocity distribution, the pattern of annual modu- 
lation of the scattering rate induced by the velocity of the 
earth can be seen as a deformed sine wave. To quantify 
it, we define the modulation as 
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where Ao and Ai stand for the amplitudes of zeroth and 
first Fourier modes of the scattering rate, respectively. 
The exact definitions of Ao and Ai and the uncertainty 
of M are studied in Appendix A. 



3 



III. COGENT DATA WITH CONSTRAINT 
FROM XENONIOO 

A. Relations between CoGeNT and XENONIOO 

The differential scattering rate measured by CoGeNT 
with 435 days data of 330 gram fiducial volume is shown 
in Fig. [U where the dots represent the spectrum with 
the L-shell electron capture (EC) background subtracted 
using the method in Ref. ^ . The errorbars include only 
the statistic uncertainties. After subtracting the L-shell 
EC background, we found that the differential rate shown 
in the bins labeled by red dots in Fig. [T] are obviously 
smaller than expected, which may be caused by some 
unknown systematic errors of the subtraction. Therefore, 
the information of these five bins is not used for fitting 
the unknown background. 
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FIG. 1: Total differential rate observed by CoGeNT with the 
prediction from the inelastic WIMP model with = 13.3 
GeV, Am^ = 33 keV, ^jv x {A^^f = 6 x IQ-^* cm^. 

As was noticed in Refs. 0, one can do a halo in- 
dependent comparison between different direct detection 
experiments. From Eqs. and (|H), the differential scat- 
tering rate can be written as 
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In the case of small nuclear recoil where we are interested 
in, since the target nucleus can be seen as rigidity, we 
have i^(|q|) ~ 1- Furthermore, if these nuclear recoils 
are generated by scattering of WIMPs, the mass of these 
WIMPs must be around 10 GeV, which is much smaller 
than the masses of germanium and xenon. Therefore, the 
reduced mass fin a in Eq. ([T^ can be replaced by M^- 
With these simplifications, Eq. (IT2|) can be written as 
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From Eq. ([3]), with the approximation ^j^j^ « M_4, one 
can see that Wmin depends only on the product ErAlA- 



Therefore, we can get 
dR pd^n ^cff 
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where one can see that the only dependence of 
dR/d{MAEr) on the target nucleus is in the factor 
A^ff/AlA- Therefore, Eq. (fT4l) can be used to translate 
differential scattering spectrum between different direct 
detection experiments. As a result, the relation between 
differential scattering rates of CoGeNT and XENONIOO 
can be written as 
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Therefore, one can get that 
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B. CoGeNT annual modulation 

The annual modulation spectra are shown in Fig. [2] 
for recoil energy in three different energy regions. One 
can see that the amplitude of modulation is still around 
15% in the region of 1.6 keVee to 2.8 keVee. Therefore, 
if we assume the only source of annual modulation ob- 
served by CoGeNT is the scattering of WIMPs, part of 
the measured spectrum in the region of 1.6 keVee to 2.8 
keVee shown in Fig. [T] should include WIMP signal. The 
quenching factor in this region is around 0.25 and the 
cut efficiency is about 0.9. From Fig. [T] we can get that 
the measured spectrum in this region is also around a 
constant which is about 3 /keVee/kg/day. Therefore, 
using Eq. ([HI) , we can get that the corresponding range 
for xenon detector is about 3.7 keVnr to 6.0 keVnr with 
differential scattering spectrum 
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Therefore, with the Poisson smearing effect [ll| , the dis- 
tribution of SI signal in XENONIOO detector can be writ- 
ten as 



Nsiin) = r]{n) 



dR 

dErPoissonin, fisi{Er))- — 
dEr 



(19) 

where rj is the acceptance of SI signal in XENONIOO, 
Psi = Ej^j-LyLcsSni/ Sec is the expectation value of the 
number of photoelectrons induced by nuclear recoil Er- 
The values of the quantum yielding efficiency Ly, the 
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scintillation efficiency L^s and the electric scintillation 
quenching factors 5ee and S'm- can be got from Ref. Q . 

In XENONIOO, the energy window of WIMP search 
region is chosen between 4-30 photoelectrons. From 
Eq. ([T9|. we can get that the spectrum in the region 
of 1.5 keVee to 2.5 keVee observed by CoGeNT is corre- 
sponding to 
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If the interaction between WIMP and nucleons is isospin- 
violating, (A^^/A^^)'^ can be as small as 1/20. There- 
fore, we can get that, if CoGeNT spectrum in the region 
of 1.5 keVee to 2.5 keVee were all induced by WIMP, 
XENONIOO should have already observed around 10 
events. However, XENONIOO reported only 3 events 
passed all cuts and only one of them lies in the low en- 
ergy region with six photoelectrons observed. Therefore, 
we can get the conclusion that WIMP signals amounts 
at most to 10% of scatterings in the region of 1.6 keVee 
to 2.8 keVee. 

However, as shown in Fig. [IJc), the modulation of the 
data is around 15%, which means either the background 
itself has an annual modulation with the same phase as 
the one generated by WIMP, or the modulation of WIMP 
is around 100% with the background having no modula- 
tion at all. In this work we pursue the second possibility 
that we assume the background has no dependence on 
time and all the annual modulation of the CoGeNT data 
is generated from WIMP. To generate 100% annual mod- 
ulation, we need to resort to a DM model in which the 
interaction between WIMP and nucleus is inelastic and 
isospin- violating as well. 

Since only a small fraction of the differential rate shown 
in Fig. [1] can be generated by WIMP. We only use our 
model to fit the modulation and use the total rate as a 
constraint. 



IV. A SIMPLE MODEL FOR 
ISOSPIN- VIOLATING INELASTIC WIMP 

Theoretically, it is not difficult to get an inelastic 
WIMP with isospin-violating interaction. The isospin- 
violating Z' model has already been proposed in the liter- 
ature [16[ . In this work, we extend this model to inelastic 
WIMP. The Lagrangian can be written as the following. 



(21) 



where (/„ and gd are couplings between Z' to up-type 
and down- type quarks, respectively, x is the candidate 
of DM stabilized by a global U{1) symmetry with mass 
m^. Now, consider the case that this C/(l) symmetry is 
softly broken by the following mass terms 
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FIG. 2: Fitting the annual modulation signal detected by 
CoGeNT. The data points represent the observed single after 
subtracting the L-shell EC background. The black curves 
show the prediction of the Z' model with = 13.3 GeV, 

Am^ = 33 keV, ctat x {A'^^f = 6 x 10"^** cm^. (a) is for 
ionization energy from 0.5 to 0.8 keVee, (b) from 0.8 to 1.5 
keVee and (c) from 1.6 to 2.8 keVee. 
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S'[x' + X' 



(22) 



where 5^ is assumed to be real. It is technically natu- 
ral to assume S'^ <C m^, since it is a parameter which 
softly breaks a U{1) global symmetry. Then, a mass gap 
Am^ = S'^ /m^ is generated between the real and imagi- 
nary parts of X- Therefore, for a 10 GeV dark matter to 
get a mass gap around 30 keV, the symmetry breaking 
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parameter S is around 20 MeV. In this work, we use xi 
and X2 to label DM and its excited state, respectively. 

In the Lagrangian shown in Eq. (f2T|) . to avoid flavor 
changing neutral current, the couplings gu and gd should 
be universal for all the three families of quarks. Then, it 
is easy to see that 



gu + 2gri 

Ap 2gu + Qd 



(23) 



where An and Ap are scattering matrix elements of 
WIMP to neutron and proton, respectively. One can 
see that if gu/gd ~ —1.125, An/Ap approaches to —0.7, 
where the sensitivity of xenon is significantly reduced 
compared to germanium. The cross section between 
WIMP and nucleons can be written as 



~p(n) 
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(24) 



It is easy to see that this U{1) symmetry is not anoma- 
lous free, however, one can always introduce heavy spec- 
tating fermions into the model to free the anomaly, and 
the upper bound of mass of the spectators is around 
64:TT^Mz' /9z' which can be far away from the reach 
of detectors. 

From the discussion of Sec. IIB, we have already seen 
that if one wants to interpret the annual modulation of 
CoGeNT as signal from WIMP, a large amount of events 
should be unknown background other than the known L- 
shell EC background. Therefore, we do not attempt to 
use the model to fit the differential rate in Fig. [1] Our 
strategy to calculate the annual modulation is the follow- 
ing. We first get an interpolating function of differential 
rate with the L-shell EC background subtracted. The 
model predicts a differential rate for CoGeNT, taking 
the green curve in Fig. [T] as an example. Then, we sub- 
tract the predicted differential rate from the interpolating 
function to get the unknown background. Assuming the 
unknown background has no annual modulation we can 
calculate the modulation predicted by the model. 

To fit the observed modulation of CoGeNT, we need 
the modulation from WIMP to be around 100%. In real- 



ity, the parameters we use are = 13 GeV, Am^ 



33 



keV, ajv x {A^^) = 6 x lO'^s cm^. The predicted differ- 
ential rate is shown by the green curve in Fig. [TJ where 
one can see that it does not exceed the observed one. 

The predicted curves of modulation for different en- 
ergy regions of 0.5 — 0.8 keVee, 0.8 — 1.5 keVee, and 
1.6 — 2.8 keVee are shown in Fig. [Altogether with the 
binned CoGeNT data. The change of the modulation 
with the energy is shown in Fig. |3l where the predicted 
data are represented with red triangles while the experi- 
mental data are shown with blue round dots. The error 
bars shown in the plot include only the statistic uncer- 
tainties, which is discussed in Appendix A. The energy 
ranges for the seven points shown in the plots are 0.5 
- 0.8 keVee, 0.8 - 1.5 keVee, 1.0 - 2.2 keVee, 1.2 - 
2.4 keVee, 1.4 - 2.6 keVee, 1.6 - 2.8 keVee, and 1.8 - 



3.0 keVee, respectively. Since the energy regions overlap 
with each other, the last six points are not independent so 
that one cannot do a analysis using the plot in Fig. [31 
The reason for choosing overlapped regions is that, the 
events are not sufficient enough. The energy bins we are 
choosing are to keep the average number of events per 
month to be not smaller than 30, so that from Eq. (jA4p 
the statistic uncertainty is smaller than 7.5%. 
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FIG. 3: Annual modulation as a function of ionization en- 
ergy. The blue dots show the modulation extracted from 
CoGeNT experiment after subtracting the L-shell EC back- 
ground, while the red curve shows the corresponding theoret- 
ical predictions from the model discussed in Sec. IV with the 
parameters chosen as the same as in the caption of Fig. [21 
The energy region of the first plot is from 0.5 to 0.8 keVee, 
the second from 0.8 to 1.5 keVee, the third from 1.0 to 2.2 
keVee, the fourth from 1.2 to 2.4 keVee, the fifth from 1.4 to 
2.6 keVee, the sixth from 1.6 to 2.8 keVee and the seventh 
from 1.8 to 3.0 keVee. The dashed green curve shows the 
corresponding predictions using the parameters from Ref 

One can see that the theoretical prediction of the an- 
nual modulation follows exactly the trend of the experi- 
mental data. The modulation is small at the first point. 
The reason is that the inelastic scattering rate diminishes 
at very small recoil energy as shown in Fig. [1] In the en- 
ergy region of 1.0 — 1.5 keVee, the predicted rate reaches 
its maximum. As a result the modulations also become 
large. In the energy region of around 1.5 keVee or larger, 
although the differential rate gets suppressed. The mod- 
ulation signal is still strong due to that higher velocity 
is needed to generate large nuclear recoil. The dashed 
green curve stands for the corresponding prediction with 
parameters used in Ref. Q , where we can see that in high 
energy region it is lower than the experimental value. 

With these parameters, the number of events that 
XENONIOO should have seen in the region of 4 to 12 
photoeletrons is about 



50 X 



/ ^off 



(25) 



In this energy region, only one event passed all the 
cuts with the expected back grou nd to be 0.6. Using 
the Feldman-Cousins method UM, to be consistent with 
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the XENONIOO result up to 90% C.L., the predicted 
number of events should be smaller than 3.6. There- 
fore, XENONIOO constraint can be avoided by choosing 
{A^i/A'^if < 0.072, which requires 

- 0.706 < An/Ap < -0.674 . (26) 

CDMS also uses germanium as target materials, which 
may render a strong constraint on CoGeNT result. Con- 
flicts from the scattering rate can no longer be alleviated 
using the idea of isospin-violating interactions. However, 
there might be issues related to the calibration of CDMS 
energy scale of the CDMS detector [19j, and the differ- 
ent temperatures in CDMS and CoGeNT may also cause 
a discrepancy between their detecting ability. Therefore, 
in this work, we are not going to take the constraint from 
CDMS into consideration. 



colliders, the upper bound on Mz' is around twice of m^. 
This upper bound is quite common in the case of isospin- 
violating, inelastic models. The reason is that, to get a 
relative enhancement of germanium over xenon, A'^^ is 
as small as 5 x 10"'^. Furthermore, in the case of inelastic 
scattering, only a very small fraction of DM particles has 
enough kinetic energy to interact with the target nucleus. 
This fraction is around 0.01 in the case with parameters 
in this work and around 0.1 in Ref. Therefore, the 
constraint from collider physics is enhanced by a factor 
of lO'* ^ 10^. Therefore, in this case, it is a generic fea- 
ture that Z' should be produced in the collider so that 
the cross section for the process in Eq. (P7)l is suppressed 
by the phase space of three body final state. 



V. COLLIDER CONSTRAINT 

One can use the process 

pp -> monojet-|-xiX2 , 



(27) 



in Tevatron to constrain this model. X2 may or may 
not decay inside the collider. However, since the mass 
gap between xi f^^d X2 is ~ 30 keV, X2 can decay to 
Xi by emitting only photons or neutrinos or some un- 
known light particles in a hidden sector. Both neutrinos 
and hidden light particles cannot be seen by the detector 
and the photons are also too soft to be detected. There- 
fore, the signal of this process is a single jet plus missing 
transverse energy (MET). 

The CDF group used data with a luminosity of 1 fb^^ 
to study this process [2l|. With the cut that both the 
transverse momentum (Pt) of the leading jet and the 
MET should be larger than 80 GeV, Pt of the second 
hardest jet should be smaller than 30 GeV and vetoing 
any third jet with Pt larger than 20 GeV, 8449 events 
were found with the expected SM background to be 
8663 ± 332. To set a 2a limit on new physics, we require 
that the cross section for the process in Eq. (|27p should 
be smaller than 0.664 pb. In this work, we work on the 
parton level process. Pioneered works have already been 
done in the literature 22| . A detailed study of Z' model 
will be published elsewhere Here, we only present 

the result of the isospin- violation case. To satisfy the con- 
dition in Eq. P^ . we can get —1.139 < gu/gd < —1.123. 
Therefore, in the study of collider constraint, we fix this 
ratio to be —1.125 which gives the strongest isospin- 
violating effect of xenon and germanium. The upper 
bounds for g^ with different values oi go are shown in 
Fig. 0] (a), with the corresponding upper bounds on 
shown in Fig. 0] (b). The sharp cliffs at around 2771,^ in 
(a) are caused by the transition of Z' from off-shell to 
on-shell. The black dashed curve in (b) shows the aN 
required to generate the modulation signal. One can see 
that to satisfy the constraints from direct detection and 
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FIG. 4: Collider constraints on isospin violating Z' model, 
(a) shows the constraint on g„ with go fix to be 0.5 (red), 
1.0 (green) and 2.0 (blue), (b) shows the corresponding con- 
straint on (Tjv. The thin dashed curves and dotted dashed 
curves are for = 7.5 and 20 GeV, respectively. The thick 
dashed black horizontal line is the value used to fit CoGeNT 
modulation in this work. The green thin dashed horizontal 
line shows ajv used in Ref. 01 . 
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VI. SUMMARY AND DISCUSSIONS 

In this work, we use a simple Z' model with inelastic 
WIMP to fit the modulation data observed by CoGeNT. 
We found that with a relatively small Wq, a 13.3 GeV 
inelastic WIMP with 33 keV mass gap from its excited 
state can be used to fit the modulation data and avoid 
the constraint from XENONIOO. We considered the con- 
straint from Tevatron on this model, and find that Mz' 
should be smaller than ^vn^. 

As noticed in Ref. [13], if Xi can only decay into w 
through kinetic mixing between Z' and S, where B is 
the gauge boson corresponding to C/(l)y gauge group in 
SM, the life time of xi '^iU be longer than the age of the 
universe and the "down-scattering" of xi on the nuclei 
will give the dominant contribution to nuclear recoil. As 
a result, the nuclear recoil tends to be larger compared 
to the "up-scattering" case so that the constraint from 
XENON 100 would be more stringent. Furthermore, in 
the "down-scattering" case, since the internal energy re- 
leased from X2 can be converted into the kinetic energy of 
the recoiled nucleus, the dependence of the nuclear recoil 
on the velocity of DM is not as strong as in the case of 
"up-scattering" . Therefore, it is not easy for the "down- 
scattering" collision to generate enough annual modula- 
tion to fit the CoGeNT data. As a result, one needs to 
invoke other mechanisms to make the decay of xi more 
rapid. 
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Appendix A: Statistic uncertainty of modulation 

The Fourier transformation of signal can be written as 

^0 = ^ ^ F(t)dt ■ 

" 1/ ^^(t)sin(27rVr)di ; 
2 

A\ = - \ F(t) cos{2T:t /T)dt ; 
^ Jo 

A, = ^Af + Af , (Al) 



where F is the observed differential rate and T is one 
year. In the analysis, to get a good statistics, each point 
in Fig. [2] represents a collection of the data in a month. 
Therefore, the Fourier analysis in Eq. (jAip should be 
discretized that 



1 ^ 

i=l 

2 ^ 

i=l 

2 ^ 

= -5]F,cos(27rzT/Ar) ; 

i=l 

Al = ^Af + Af , (A2) 



where N is the number of partition in one year. 

If all the signals observed are induced by noise, 
one may still observe an annual modulation using the 
Eq. (fTT|) . In this case, we can assume Fi in Eq. (jA2p 
obeys a Poisson distribution with expectation value F . 
Then, it is easy to get that 



(A^) = {AD = , 
(Af ) = (Af ) = 2F/N . (A3) 



Therefore, the statistical uncertainty of the modulation 
can be written as 



5M ^ \ ^ . 
V NF 
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